SELEN: A Fortran 90 program for solving the "sea-level equation"

We present SELEN, a Fortran 90 computer program for solving the “sea–level equation”, the law that governs the sea–level variations associated with the melting of the Pleistocene ice–sheets. As shown here by a number of examples, SELEN also allows to study a suite of geophysical processes accompanying the mass redistributions associated with glacial–isostatic adjustment, such as vertical deformations of the solid surface of the Earth, and variations of the shape of the geoid. SELEN is built in a simple and easily testable manner, and requires modest computer resources. The source code, which is freely available on request, can be easily modified to implement new features of the sea–level equation or to focus on specific geodynamical aspects of glacial–isostatic adjustment.